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, We show in this paper that the observation of the angular distribution of upward-going muons and 

^SJ ' cascade events induced by atmospheric neutrinos at the TeV energy scale, which can be performed 

, ^ ' by a kilometer-scale neutrino telescope, such as the IceCube detector, can be used to probe a large 

^ ^ . neutrino mass splitting, |Am'^| ~ (0.5 — 2.0) eV^, implied by the LSND experiment and discriminate 

among four neutrino mass schemes. This is due to the fact that such a large mass scale can promote 
non negligible —> i/e, z^t/S'm — *■ f'e, i't conversions at these energies by the MSW effect as well as 
lO vacuum oscillation, unlike what is expected if all the neutrino mass splittings are small. 

(N : 

PACS numbers: 14.60.St,14.60.Pq,14.60.Lm,95.55.Vj,95.85.Ry. 
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O ■ I. INTRODUCTION 

(N : 

■ Atmospheric and solar neutrino experiments present today compelling evidence of neutrino flavor oscillations in- 
' duced, respectively, by the mass squared differences |Ato^.j,„| ^ 3 x 10""^ eV^ and Atoq ^ 7 x 10~^ eV^ 2]. These 
, results have been independently supported by the K2K experiment Q, confirming the v^/i>^^ conversion observed in 
atmospheric neutrino data, and by the KamLAND experiment Q, recently observing Df. disappearance compatible 
with what has been seen by the solar neutrino experiments. 

In the standard three neutrino flavor oscillation framework, the above results imply neutrino oscillation lengths 
Q \ at energies above about 0.3 TeV much greater than the Earth's diameter, strongly suppressing Ve/^r and 

1 v^jv^ oscillations. Moreover, the ratio between v^jVf^ and Vp^jv^ atmospheric fluxes drops as their energy 

^ ' increases so that atmospheric neutrino data above 0.5 TeV are expected to be basically composed of v^^ and v^, ja]. 
• ^ : The so called prompt atmospheric neutrinos coming from the decay of charmed particles as well as extra-galactic 
v^, Vf./vr, Vt contributions are very small compared to the conventional atmospheric v^lVfi fluxes (which mainly come 
^ from the decay of charged pions and kaons) in the 0.5-10 TeV energy range p|. 



However, as we will see, atmospheric v^, Pe/j^T, Vr fluxes in the TeV energy range can be non-negligible compared to 
the Vfj,/Vfj, if the result of LSND experiment \j\ which indicates — > v,. transition generated by neutrino oscillations 
with a large mass scale, jAm^g^^l ^ (0.5 — 2.0) eV^, turns out to be correct. The LSND result, which will be very 
soon tested by the MiniBooNE experiment @, if conflrmed, can be interpreted in combination with the atmospheric 
and solar neutrino data as a signature of an additional neutrino, Vs, sterile in nature. The minimum scheme necessary 
to explain at same time all three indications of neutrino oscillations involves four neutrino mass eigenstates . We 
are aware that current analyses performed with the entire set of neutrino data from accelerators, reactors, solar and 
atmospheric neutrinos seem to disfavor four neutrino mass scenarios See, however, Refs. |0 which casted 

doubt on conclusions drawn in Refs. lll|. More importantly, no experiment until now has conflrmed or ruled out 
the LSND result and one should wait for results from MiniBooNE. 

We will show in this letter that in the framework of four neutrino mass schemes, atmospheric neutrino data at the 
TeV energy scale may contain a non negligible amount of Ve or Pe, due to the interplay of the Earth matter with the 
LSND mass scale through the so called MSW effect , and of Vr and Vr , due to almost vacuum oscillations induced 
by this same mass scale. To our knowledge it has been never mentioned before in the literature that the LSND mass 
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FIG. 1: Four neutrino mass schemes used in this paper. 



splitting could give rise to resonant or — *■ D,, conversion for atmospheric neutrinos. These I'e, Pe/'^r, i't 

events can be observed by a large neutrino telescope, such as the proposed IceCube detector by measuring the 
angular distribution of upward-going muons and cascade events. This is because the zenithal dependence of upward- 
going muons as well as of cascade events will be modified by these new contributions with respect to what would be 
expected for the standard three neutrino oscillation framework. The zenithal distortions can not only discriminate 
among the different types of four neutrino mass schemes, but also they are sensitive to the sign of the LSND mass 
splitting, opening a new window to test four neutrino mass scenarios. 



II. FRAMEWORK OF FOUR NEUTRINO MIXING SCHEMES 

We will work in the minimal model structure needed to describe the three evidence for neutrino oscillations (atmo- 
spheric, solar and LSND): a four neutrino mass scheme 0|. We will adopt the same parametrization and some useful 
notations used in Refs For what concerns us here, the sign of Am| and of Ato^^^, are irrelevant, only the sign 

of AmLSND will be significant. We can then classify the possible four neutrino mass schemes into two categories: (a) 
(3+1), when a single neutrino mass state is separated from the other neutrino mass states by a large mass splitting 
(IAtolsndI); and (b) (2+2), when two pairs of neutrino mass states are separated by a gap given by the lAmLg-^ol 
scale. Moreover, when 1/4 is a heavy (light) state we will call it normal (inverted) mass hierarchy. We show in Fig. ^ 
the (2+2) and (3+1) mass schemes for both, normal and inverted, mass hierarchies which will be considered in this 
paper. 

The evolution equations for neutrinos in the fiavor basis are 



"f. (1) 



where z// = [veVf^Vrfsy and is the neutrino energy. We define I^m1^^„ = Am^j and Am^^^, — Am^j, with 
Amfj — mf — 771^, for i, j = 1,2, 3, 4. The diagonalized neutrino mass squared matrix is 

= diag(0, -Aml,^„ + Am|, Am^^,, + Am| - Am^,^^, Am|), (2) 

with Am| — Amli for {2 + 2) cases and 

M2 = diag(0, Am2 , Ami + Am^^,,, Am^,^^ + Am^ ), (3) 

with Am| = Am2i for [3 + 1) ones. We use for the 4x4 mixing matrix U the same parametrization as in Ref. (15i |. 

U ~ R34R24R23R14R13R12, (4) 

where Rij are rotation matrices in the subspace (i, j) with angle 9ij. The matter potential matrix is 

A = diag(Acc,0,0,Anc), (5) 

with Acc — and A^^ = ^cc(l ~ ^e)/(2ye), where Gp is the Fermi constant, Yg the electron fraction and p 

the matter density. For antineutrinos evolution equations similar to Eq. Q exist but with A — s- ~A. We have ignored 
all CP violation phases, as they do not play any important role in our study. 
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In order to make clear the connection among mixing parameters and experimental observations, we will introduce 
some useful definitions. Let us first define the quantities da {a — e, /i, r, s) as 



and rja (a — e, fi, t, s) as 



da = 



|J7q4P for (3+1) mass scheme, 

\Uai? + \UaA? for ('g+g) mass scheme, ^ ' 



|C^QiP + |C^a2p for (^5+ mass scheme, 

|C^Qi|^ + |C^a4|^ for mass scheme, ^' 



which can be related to the neutrino survival and oscillation probabilities at short and long-baseline experiments as 
we will see below. 

Note that r^s measures the sterile admixture in solar neutrino oscillations, this parameter describes the fraction of 
sterile neutrino participating in solar oscillations: 7]s — Q (1) corresponds to pure active (sterile) Ue ^ (^^e ^ i^s) 
oscillations. Solar neutrino data prefer a small admixture of sterile neutrino in the Attiq scale, i.e., ijs <C 1. The 
sterile admixture for atmospheric neutrinos is characterized by the parameter ds, in the Am^rj,,^ mass scale, dg = 1 (0) 
corresponds to pure active (sterile) Vr {Vfi — > Vs) oscillations. Atmospheric neutrino data prefer a small 

admixture of sterile neutrino in the Am^^^, scale, i.e., 1 — <C 1. 

We can relate da to the neutrino survival probability in vacuum for short-baseline disappearance experiments, 
ignoring ©(Am^^^j) corrections, by 

=1-A„,asin2(^^^i^ =l-4d„(l-d„)sin2(^^^i), (8) 

where L is the baseline distance. From the null results of short-baseline experiments Ae;e and A^.^^^ are limited to be 
very small Q. 

The conversion probability, relevant for short-baseline appearance experiments, can be written as 

P.^^.^=A,,, sin^l^^^i), (9) 

where 

4|l7e4p|/7p4p for f5+ J j mass scheme, , ^ 

4|f7eiC/;i + UedJ*^4,? for (2+2) mass scheme. ^ ' 



Here the amplitude A^^je can be identified as the LSND mixing amplitude sin^ 20lsnd 0- 

For the long-baseline reactor (L < 1 km) and accelerator experiments [L < 1000 km), we have to include effects of 
the Am^TM niass scale, and the conversion probability, P™.^^^ becomes 

P-_.,^ = 1 _ 2 d.(l - 4) - sin^ (^^i) , (11) 



with 

Ac 



4rye(l - de — rje) for (3+1) mass scheme, 

4 I ?7e2 P I C^e3 P for ( 2+2) mass scheme. 



(12) 



To accommodate the {2+2) mass scheme, d^, Acn and A^.^e have to be small and d^ close to 1 (to keep A^^ 
amplitude small). Also the parameters rjs and dg are not independent, rjs = dg as noted by the authors of Ref. |l6| . 
From Ref. |llfi| we know that the global fit of atmospheric and solar neutrino data prefers either the limit r}s = d^ ^ \ 
or the limit rjs ^ ds <ti 1. For the {3+1) mass scheme, d^,de, Ach, ^/^;e, ris and 1 — d^, all, have to be small 

In this paper, in order to give an explicit example of the effect of four neutrino mass schemes we will fix the values of 
the mixing parameters. These were taken within the parameter region allowed by solar and atmospheric neutrino data 
as well as by the null results of reactors and accelerators (with the exceptions of K2K 3] and KamLAND Q) 0,1^- 
We have scanned a grid of values to find the maximal possible values of the six mixing angles in order to maximize 
the oscillation effect driven by the LSND mass scale in IceCube, using the constraints on da, A^-e, Acu, rje and rjs 
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FIG. 2: Neutrino and anti-neutrino probabilities, for diametral trajectories {cos6z = —1), in the (2+2) scheme with normal 
hierarchy for Amlsm = 0.5 eV^ without (left panels) and with (right panels) matter effects taken into account. We have used 
AmATM = 3 X 10~^ eV^, Attiq = 7 x 10~^eV^ and the mixing parameter as given in Table|l] The standard three generation 
oscillation probabilities are labeled as 3(g). 

shown in Ref. [l^. The result of our scan is shown in Tabled For the neutrino mass spUttings, we have used the best 



3 X 10 eV^ and Am?, = 7 x 10 ^ eV^. For definiteness we also have set lAr 



= 0.5 eV^ 



We have numerically solved Eq. to compute P(j^^ Va) for a = e, /i,T and s, the oscillation probabilities that 
are relevant for atmospheric neutrino data. Let us make some observations that can help the reader to understand 
our numerical results. 

(2+2) case: 

From the expressions of the evolution equations, one can expect matter enhanced MSW resonant conversion in 
Earth for the neutrino (antineutrino) channel u^^ — > {Vfi I'e), for the 2+2 inverted (normal) hierarchy if the 
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Case 


S\V? 9l2 


• 2 a 

Sm 0l3 


sin^ 6114 


sin^ 6*34 


sin^ 024 


sin^ 623 


77s 


ds 


2+2 


0.001 


0.01 


0.26 


0.8 


0.03 


0.5 


0.2 


0.2 


3+1 


0.26 


0.035 


0.02 


0.45 


0.045 


0.5 


0.13 


0.52 



TABLE I: Values of the mixing parameters used in this paper. 



following condition is satisfied, 

^ |A<.„I X . TeV X ,4 - X (^) X (aM£f ) . ,,3, 

Other resonances can be induced by the Am^^^, mass scale at E^, < 100 GcV, but they are irrelevant for our purpose 
here. In the (2+2) mass scheme, de ~ 1 and ^ 1 so that the resonant condition can be achieved in the TeV energy 
range. For the — s- v-^ channel we expect almost vacuum oscillations. 

In Fig. 12 we show, for diametral trajectories (cos^z = — 1), the behavior of the neutrino and antineutrino v^j, — > 
and v-TTVei transition probabilities for the (2+2) case with normal hierarchy. To provide a better understand- 

ing of what is going on, in panels (a)-(d) the matter effects were not taken into account while in panels (e)-(h) they 
were. We also show for comparison what would be expected for the standard three flavor neutrino oscillation scenario 
(labeled as 3(g)). 

If we compare Fig.|2Ic) and Fig.|2Ig), we see a sizable De conversion probability when matter effects are taken 

into account. This resonance behavior is due to the LSND mass scale and Earth matter effects Looking closely 
to Fig.|2Ig), we can see three peaks: one at ~ 0.57 TeV for a typical value of Earth core density (pcore ~ 12g/cc), 
and another for E^ ~ 1.5 TeV for a typical value of Earth mantle density, (jOmantic ~ 4g/cc). A third peak appears 
for E^ - 0.8 TeV and its origin is due to a parametric resonance. Similar behavior was noticed sometime ago in 
the context of neutrino oscillations induced by Am^^^, for z/^ Vs To determine the behavior of oscillations in 

— > Vr channel, we can use Eq. lO, changing e — > r, so that P^^^^y^ — Af^-T sin^ {Aml^j^^L / (iE,^)) . The conversion 
amplitude for this channel, A^;r = 4 cos^ ^^34 cos^ 6*24 sin^ ^24 depends only on sin^ ^^24 and sin^ ^^34. We confirmed 
numerically that this is correct and we have virtually vacuum oscillations induced by Ato^snd as you can see comparing 
Figs.|2Ib) and (f). In the case of inverted mass hierarchy, the neutrino and antineutrino probabilities are interchanged. 

(3+1) case: 

In the (3+1) mass scheme, we cannot expect a strong resonant Ve {v^ — > Ve) transition due to the LSND 

mass scale at TeV energies. This can be explained as follows. First, both and {Pf_i and D^) are basically 
distributed among the mass triplet so that we cannot expect a significant direct transition between them driven by 
Amlgy^j^. However, a second order transition can be achieved via which is mainly composed of z/^ as Vfj, I's v^, 

{^t^ ^ Vs ^ Ve)- 

For the normal mass hierarchy, the first step, — s- {Vfj, Pg), occurs via almost vacuum oscillation (a MSW 
resonance effect), while the second step, Us — » Ve {vs i^e), occurs via a MSW resonance effect with E^ given by 
Eq. (|13|l with de — replaced by ds (almost vacuum oscillation) . The strength of these transitions is regulated by 
the value of On and 024, which are very much constrained by current data (see Table I). This is why we do not expect 
to observe a large number of (Pe) events in the (3+1) scheme. In the case of inverted mass hierarchy, the neutrino 
and antineutrino probabilities are interchanged. 

On the other hand, (P^ — > transitions driven by the LSND mass scale can occur directly mainly via 

vacuum-like oscillation, since there is a significant amount of Vr in the v^, state. This is due to the large value of ^34 
(see Table I), which is the maximal value still allowed by data. 

These observations are confirmed by our numerical calculations. In Fig. |31 we show, for diametral trajectories 
(cos^z = — 1), the behavior of the neutrino and antineutrino v^^ — )■ i/^ and i/^ ~^ v-T^Vf^^Vg transition probabilities for 
the (3+1) case with normal hierarchy in a similar way as in Fig. |2] 

From these observations we should expect at the TeV energy scale Vf, matter enhanced resonant transitions 

for neutrinos (antineutrinos) if the LSND mass hierarchy is normal (inverted) in the (2+2) scheme and — s- as 
well as — > Vr vacuum-like transitions in both mass schemes. 
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FIG. 3: Neutrino and anti-neutrino probabilities, for diametral trajectories {cos6z = —1), in the (3+1) scheme with normal 
hierarchy for AmtsND = 0.5 eV^ without (left panels) and with (right panels) matter effects taken into account. We have used 
'^^ATM = 3 X 10~^ eV^, Am0 = 7 x 10~^eV^ and the mixing parameter as given in Table|l] The standard three generation 
oscillation probabilities are labeled as 3(g). 



III. TYPES OF SIGNAL EVENTS 



We have studied the observable consequences of these four neutrino schemes at high energies, in a future neutrino 
telescope. For definiteness we have assumed a kilometer-scale neutrino telescope as the one proposed by the IceCube 
collaboration 14] that, for simplicity, has been simulated with perfect efficiency as well as perfect angular and energy 
resolution. We have supposed that the IceCube detector can be described by a cylinder with 900 m of height and 
500 m radius. 

We have calculated two types of measurable signals in this detector: upward-going muons and cascades events. The 
calculations were done using the atmospheric neutrino fluxes, denoted by <i>^^ , a = e, /i, of the Bartol group ^'^d 
the cross sections from Ref. We have also included in our calculations the effect of neutrino absorption when 

crossing the Earth, taking from Ref. [T^. 
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A. Upward-going muons 



Muon neutrino charged current interactions, occurring inside the Earth close enough to the detector, may produce 
energetic muon events which can cross the detector volume, giving rise to the so-called upward-going muons. These 
events are characterized by a very good angular resolution. 

We have calculated the number of upward-going muon events as a function of the zenith angle, Nf^(6z), as in 
Refs. mils, 

N^{0.) r e,)'^^Al^^^dE^ + antineutrinos, (14) 

•^Bumin atj^^a COS 



with 



dEf^d cos dz 



= Na 



poo pOO pOO 

/ dEfj^o / dE^ / dh K„^{h,cos9z,E^) 
Je^ Je^o Jo 



d^,^{E,,ez) da{E,,E^o) p wrp n\ (^^^ 

X , ^Py^f^ ^ ^m) TT; i'rock[Ef,o,E^)b{E^,Oz), (15) 

dE^dcosOz dEfj^o 

where T is the livetime that we assumed to be 1 yr, Na is the Avogadro number, i?^o is the energy of the muon produced 
in the neutrino interaction and E^ is the muon energy when entering the detector, Frock{E^o, i?^i) is a function that 
takes into account the muon energy loss from the point where the muon is created to where it is detected |20|, Ka 
is the slant distance distribution normalized to one p^ . Pli^f^ — > f^j) is the muon survival probability, A(Rmin, Oz) 
is the effective area for a muon traveling a distance i?min inside the detector coming from the direction 9^ |20j, a is 
the i'^,N fiX, neutrino cross section, cosOz labels both the neutrino and the muon directions which to a very good 
approximation at the relevant energies are coUinear, S{E^,62) is a shadow factor that is related to the attenuation 
of neutrinos traversing the Earth [l9l |. For a given zenith angle, the threshold energy cut is obtained by equating 
R{Efj,.min) = Rmin — 100 m, where i?(-E'^,min) IS the muon range function. 

In Fig. 2] we show the zenithal angle distribution of these events expected in the case of the standard three flavor 
oscillations (or no oscillations at all) and for each of the four neutrino mass schemes assuming |Am^g„n| — 0.5 eV^ and 
the oscillation parameters as given in Table ^ It is clear from this plot that the distribution of upward-going muons 
reflects the v^/v^^ disappearance into other flavors promoted by each of the four neutrino schemes, making them quite 
distinguishable from the standard three flavor oscillation case. It may be even possible to distinguish among (2+2) 
and (3+1) schemes, although their normal and inverted cases are virtually the same due to the fact that the mean 
energy of these events is much lower than the resonant energy [E^ ~ 0.8 TeV for /^ml^^^ = 0.5 eV^). From the 
Eq. (fnji we note that all muons with an energy sufficient high to travel a distance i?min are classified as upward-going 
muons, the original muon energy is integrated out so that the correlation between the original neutrino energy and 
the final muon energy is weak. 

As a test of capability to discriminate the four neutrino mass schemes from the standard three generation oscillation 
we performed a simple analysis, as follows. If we define N^{Ai/) (A'^(3g)) as the number of upward-going muons 
in the i-th zenith angle bin in four neutrino mass schemes (standard three neutrino mass scheme) showed in Fig. 0] 
then we can compute the function 



where (^[^i,) (^(3g)) statistical error on TV* (4i/) {N^{3g)), for simplicity we have ignored systematic uncertainties. 

The parameter 77 is a normalization factor. If we assumed predictions of upward going muons in standard three 
generation are the correct ones, then the only possible change is a global normalization factor. We scan the parameter 
1] to look for the best fit and for any of neutrino mass schemes the minimum is much higher then the degrees of 
freedom, 20-1 (20 points minus 1 normalization factor), with a goodness of fit well below 1%. Then we conclude that 
for the choice of mixing parameters in this paper, it is very likely that the zenith angular distribution of upward going 
muons is sufficient to be separated from the usual standard three neutrino oscillations. From this we can say that 
independent of which four neutrino mass scheme nature chooses, we can always have a distinctive zenith distortion 
in upward-going muons. 
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FIG. 4: Expected zenithal distribution of upward-going muon events in the IceCube experiment after one year of data taking 
for the standard three flavor oscillations (3g) and the different four neutrino mass schemes. 



B. Cascades 



Electromagnetic and hadronic showers can be produced by ionizing products of the neutrino interactions with the 
ice. These types of events can be created by electron neutrino charged current interactions {v^ ~ CC), tau neutrino 
charged current interactions {v-r — CC) and neutral current interactions {v — NC). These events, which are called 
cascade events, are characterized by a very good energy resolution. In principle, muon neutrinos also produce a 
hadronic shower and a muon and this case contribute to cascade events. We assume events with a muon track always 
can be separated from the other cascade events and we will not include them in our computation. We have calculated 
the number of cascade events for a given zenith angle, A^cascadc(^'z), as in Refs. [23. l2^. 

iVcascadc(^?.) = TVNa [ '^^^J^T^^IT^ P ^ ^e) V ME,)dE, + antincutrinos, (17) 

/n n h,..a rcif^H ^ '—^ 
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here V is the IceCube experiment volume (1 km'^), fi{Ei,,z), i 
showers), v—NC, are functions which depend on the interaction, 



Ve—CC, Vr—CC (hadronic plus electromagnetic 



d(7cc{Ey,y) 



^shr,inin ; , 



d(Jcc{Eu,y) dn{Er) 



Q{E,{y + (1 - y)(l - z)) ~ S.hr,mi„), 



fl^^cciEu) 



dy dz 

, I , ,dacc{E^,y) dn{Er) . . , ^ /n ^ n 



■"'-^^^^ J J - dy dz 
U^^c{E.) = Jdy^^^^^^^^eiE^y-E.. 



'shr,inin ) i 



for - cc 
for Vr — CC hadronic 
for Vr — CC electromagnetic 
for v-^C 



where Ej- is the tau energy, y, z, z' are the original neutrino energy fractions carried by the relevant lepton/hadron 
at play, and we have considered only events which have a total shower energy greater than 500 GeV (i^shr.min = 500 
GeV). 

For Ur—CC events we can have hadronic and electromagnetic cascades, respectively, from z^t — > t — *■ e decay 
channel and from Vr t h, where h means any hadrons. For v—NC {vg—CC) events we have respectively 
hadronic (electromagnetic) showers. The functions dn{Er)/dz and dn{Er) / dz' are given in Ref. |2^. Although we 
also integrate on the neutrino energy as for the upward-going muons computation, the average shower energy for this 
sample is more closely related to the average initial neutrino energy. 

Since for the energies considered here {E^ > 1 TcV), the atmospheric neutrino fluxes are basically composed of i/^ 
and and according to the usual three neutrinos oscillation scenario electron and tau neutrino appearance are 
suppressed, therefore one generally expects that only the neutral current contribution can produce cascades at these 
energies. A possible exception to this would be perhaps charged current contributions induced by Ve/i^r neutrinos 
emitted in microquasar jets psj . but these events are very well collimated with the source and can be easily separated 
from the atmospheric neutrino data. As we have shown earlier, in four neutrino schemes, electron and tau neutrino 
components will appear at these energies so that an increase on the number of cascade events should be expected. 

We show in Fig. \E\ the zenithal distribution of cascade events that can be expected in IceCube after one year of 
data taking, for the standard three neutrino oscillation case as well as for all the four neutrino schemes, assuming 
|Ato^sndI = '-'■5 eV^ and the oscillation parameters as given in Table |I| All four neutrino cases, for these values of 
the parameters, can be distinguished from the standard case. We observe that there is an increase on the number of 
cascade events at cos^z ~ —1 for the (2+2) inverted scheme due to the resonant appearance of Ve similar to what 
is shown in Fig.[5](g) for Vf.. There is no such a prominent peak for the (2+2) normal scheme due to the fact that 
in this case the resonance is in the antineutrino channel (where one has a lower flux and smaller cross section). The 
increase of the number of cascade events in the (3+1) schemes is almost exclusively due to Vr appearance. In fact, 
the angular distribution of cascade events reflect the flavor content of the atmospheric neutrino flux as measured at 
the detector site. 

We observe that cascade events are more effective in separating the mass hierarchy of four neutrino schemes than 
the upward-going muon events. This is not surprising since the 500 GeV energy cut imposed on the cascade showers 
selects events closer to the resonant energy where the effect of the large IAtolsndI = 0-5 eV^ is more sizable. For 
higher values of |AtolsndI the differences among the standard and the four neutrino oscillations distributions could be 
even larger than what we have shown in Figs.0]and|Slif not for the fact the mixing angles are much more constrained 
by data for higher masses In any case the general behavior of the curves shown in this paper will still hold for 
IAtolsndI as high as 2.0 eV^. 

We note that unlike the case of upward-going muon events, cascade events can be observed also for cosOz > 0, 
which correspond to down-going events. For such events, for any four neutrino mass scheme we have considered here, 
we do not expect any significant deviation of the zenith angle distribution from what is expected in the standard 
three neutrino oscillation scheme. In the standard scheme, in good approximation, the zenith angle distribution is 
symmetric with respect to cos dz = 0. This is because of the fact that for down-going cascade events, there is no matter 
effect and the distance traveled by neutrinos is much smaller than the typical oscillation length expected for jAm^sfj^l 
scale and E^ ^ 1 TeV so that no significant Ve, v^lv^, Vt induced events are expected. Therefore, by comparing the 
zenith distributions for cos^z > and cos^z < 0, it is expected that we can discriminate four neutrino mass schemes 
from the standard three one with better sensitivity. 

Finally, let us note that the uncertainty in the prediction of the atmospheric muon neutrino flux, which is currently 
about 15%, could potentially spoil the discriminating power among the different signals. This can be, to some extent, 
overcome by combining cascade and upward-going muon events, and so reducing the uncertainty in the overall ffux. 
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FIG. 5: Expected zenithal distribution of cascade events in the IceCube experiment after one year of data taking for the 
standard three flavor oscillations (3g) and the different four neutrino mass schemes. 



Also, the shape of the event distributions, here the key point in disentangling the different signals, is not very much 
affected by the overall flux normalization. 



IV. CONCLUSIONS 



We have studied the role that can be played by a large mass splitting implied by the LSND signal, |ATO^g„^| ~ 
(0.5 — 2.0) eV^ in a four neutrino mass scheme for atmospheric neutrino data in the energy range ^ 1 TeV. We have 
shown that upward-going muon events as well as cascade events, that can be measured by a kilometer scale detector 
such as the one proposed by the IceCube collaboration, can indicate the presence of the fourth neutrino through 
the production of Ve , and Vt , Dr which are not expected in the standard three generation oscillation case where 
we have only the small mass splittings Attiq and Am^.[,^,. Furthermore the zenithal dependence of cascade events is 
even sensitive to the sign of Atolsnd because the matter effect play an important role, leading to resonant v^. or 

1?^ — > Ve conversion depending on the sign of Ato^snd ■ We note that the standard three generation oscillation scheme 
cannot induce these signals of Vel^T appearance at the 1 TeV scale, and therefore, they constitute an unique signature 
of a large mass scale present in a four neutrino mass scheme. 

A final remark. If MiniBooNE observes a positive signal, the effects shown here should be expected independent 
of the four neutrino mass scheme. On the other hand, if MiniBooNE results turns out to be negative this may be 
viewed as a complementary test of the presence of a large neutrino mass splitting in nature. 
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